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Abstract 

Background: Prophylaxis to prevent relapses in the central nervous system after childhood acute lymphoblastic 
leukemia (ALL) used to consist of both intrathecal chemotherapy (CT) and cranial irradiation (CRT). CRT was mostly 
abolished in the eighties because of its neurotoxicity, and replaced with more intensive intrathecal CT. In this study, 
a group of survivors treated with CRT before 1983 and another group treated without CRT thereafter are 
investigated 20-25 years later, giving a much stronger perspective on long-term quality of life than previous 
studies. The outcomes will help to better understand these groups' current needs and will aid in anticipating late 
effects of prophylactic CRT that is currently applied for other diseases. This study evaluates oscillatory neuronal 
activity in these long-term survivors. Power spectrum deviations are hypothesized to correlate with cognitive 
dysfunction. 

Methods: Resting state eyes-closed magnetoencephalography (MEG) recordings were obtained from 14 ALL 
survivors treated with CT+CRT, 18 treated with CT alone and 35 controls. Relative spectral power was calculated in 
the 5, 9, al, a2, (3 and y frequency bands. The Amsterdam Neuropsychological Tasks (ANT) program was used to 
assess cognition in the executive functions domain. MEG data and ANT scores were correlated. 

Results: In the CT + CRT group, relative 0 power was slightly increased (p = 0.069) and ct2 power was significantly 
decreased (p = 0.006). The CT + CRT group performed worse on various cognitive tests. A deficiency in visuomotor 
accuracy, especially of the right hand, could be clearly associated with the deviating regional 9 and a2 powers 
(0.471 < r< 0.697). A significant association between decreased regional ct2 power and less attentional fluctuations 
was found for CT + CRT patients as well as controls (0.078 < r< 0.666). Patients treated with CT alone displayed a 
power spectrum similar to controls, except for a significantly increased level of left frontal a2 power (p = 0.030). 
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(Continued from previous page) 

Conclusions: The tendency towards global slowing of brain oscillatory activity, together with the fact that 
dementia has been reported as a late effect of CRT and the neuropsychological deficiencies currently present, 
suggest that the irradiated brain might be aging faster and could be at risk for early-onset dementia. The CT group 
showed no signs of early aging. 

Keywords: Late effects, Childhood cancer, Magnetoencephalography, Resting state, Oscillatory power, 
Neuropsychology, Accelerated aging 



Background 

Acute lymphoblastic leukemia (ALL) is the most com- 
mon malignancy diagnosed in children, representing al- 
most one third of all childhood cancers. There is a peak 
incidence at 2-5 years of age. Over the last decades, cure 
rates for ALL patients have gone from less than 5% to 
over 80% due to improved treatment protocols [1], An 
essential part of ALL treatment is central nervous system 
(CNS) prophylaxis, mostly performed by cranial irradi- 
ation (CRT) and/or intrathecal chemotherapy (CT). In 
the Netherlands, CRT was abolished in 1983 [2]. 

As survival rates for ALL continue to rise, the focus of 
clinicians is increasingly shifting to the late effects of 
therapy and related quality of life [3]. CRT patients have 
a 13.6 times higher chance of developing secondary neo- 
plasms in the irradiated field within 20 years following 
treatment [4]. Late neurotoxic effects of CRT are also 
well described. Many studies have reported a causative 
link between prophylactic CRT and long-term neurocog- 
nitive deficits among ALL survivors, for example mem- 
ory problems, cognitive slowing and attention deficits 
[5-8]. Cognitive effects beyond 10 years after treatment 
are largely unknown [9]. This study provides data from 
patients, on average 25 years post childhood treatment, 
giving a much stronger perspective on long-term quality 
of life than previous studies. These outcomes not only 
help to better understand these groups' current needs, 
but can also aid in anticipating late effects of currently 
applied prophylactic CRT that is still part of treatment 
for ALL in some countries, but also for brain tumours 
and small cell lung cancer. 

Neurocognitive late effects of chemotherapy alone are 
relatively subtle compared to the effects of CRT. How- 
ever, a significant percentage of children that received 
chemotherapy alone also suffer from neurocognitive def- 
icits, particularly in the domain of executive functioning 
(EF) [10-13]. The underlying mechanisms of neurotox- 
icity, caused by either CRT or CT, are not yet fully 
understood. Structural imaging has revealed smaller 
white matter volumes in this population, which corre- 
lated with cognitive dysfunction [9]. Functional neuroi- 
maging methods like functional magnetic resonance 
imaging (fMRI), electroencephalography (EEG) or mag- 
netoencephalography (MEG) have rarely been applied in 



late effect studies of cancer treatment, while these meth- 
ods are very suitable to study subtle brain dysfunction. 
Among the available functional imaging techniques, 
MEG has the unique property of combining high tem- 
poral resolution with good spatial resolution. The MEG 
signal is quantified in spectral power, which indicates 
the strength of magnetic induction fields generated by 
oscillatory neuronal activity and is decomposed into fre- 
quency bands [14]. Relative power is the percentage of 
power in any band compared with the sum of the power 
of all frequency bands. A power spectrum that displays a 
shift towards relatively more power in the lower fre- 
quency bands is generally considered to be a patho- 
logical sign of slowing brain activity and is associated 
with deteriorated cognition [15,16]. Also, increased 8 or 
6 power is known as a nonspecific sign of brain path- 
ology [17]. In this paper we will describe relative spectral 
power changes assessed by MEG during resting state 
with eyes closed. To our knowledge, this is the first 
study to use MEG in this population. We hypothesize 
that the pattern of oscillatory activity will deviate from 
that of controls and that the deviations will correlate 
with cognitive dysfunction. 

In 1989, Tucker et al. conducted a study of long-term 
neurocognitive effects (11.5 ±7 years after treatment) of 
CT + CRT in ALL and Non Hodgkin Lymphoma (NHL) 
survivors treated at age 15-69 [18]. They compared the 
effects between a group treated at 25 years of age or 
younger and a group treated at an older age. In the 
younger group, abnormalities were more outspoken, indi- 
cating higher susceptibility to neurotoxicity at a younger 
age. Using EEG, they found excess 8 activity and slightly 
slowed a. The abnormalities were interpreted as minor, 
but relatively more severe in the younger group, reflect- 
ing subtle cerebral cortical dysfunction. Therefore, an 
association between worse outcome and younger age at 
treatment is hypothesized. 

Methods 

Subjects 

Seventy-five eligible ALL survivors were identified from 
patient records of the VU University Medical Center 
(VUMC) and the Academic Medical Center (AMC) 
Amsterdam. They were sent a letter explaining the study 
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and were subsequently contacted by phone. Unfortu- 
nately, 45% were unwilling to participate or met one or 
more of the exclusion criteria (i.e. use of centrally acting 
drugs, active psychiatric disease or symptoms, pre- 
existing CNS disorders, metal parts in the body or preg- 
nancy). Survivors who were willing to participate were 
requested to recruit a control (sibling, partner or friend, 
n = 44). The ALL survivors were treated according to 
two different protocols developed by the Dutch Child- 
hood Leukemia Study Group (DCLSG). The first group, 
CT + CRT (n = 18) was treated according to DCLSG 
protocol ALL-5 (1979-1984), which consisted of stand- 
ard dose CT (intravenously administered vincristine, 
prednisone, L-asparaginase, daunorubicine, intrathecal 
injections of methotrexate (MTX) and prednisolone) 
and CRT (2500 Gy). The second group, CT (n = 23), was 
treated according to DCLSG protocol ALL-6 (1984- 
1988), which consisted of standard dose CT (intravenous 
vincristine, dexamethasone, L-asparaginase, intrathecal 
MTX and prednisolone), without CRT. Duration of 
treatment was approximately 2 years. 

Two CT + CRT patients were excluded from data ana- 
lysis because meningiomas were discovered during as- 
sessment. Additionally, data from two CT + CRT patients, 
five CT patients and nine controls could not be used due 
to artefacts discussed in the magnetoencephalography 
section below. None of the CT patients or controls pre- 
sented with CNS disorders. Eventually, data from 14 CT + 
CRT patients, 18 CT patients and 35 controls were used 
in the analysis for this paper. Characteristics of these par- 
ticipants are reported in Table 1. 

The ethical principles of the Helsinki Declaration were 
followed and approval was obtained from the Medical 
Ethics Committee of the VU University Medical Center 
(n° 2006/200). All participants signed an informed con- 
sent form. 

Magnetoencephalography 

MEG was recorded for ten minutes while subjects were 
inside a magnetically shielded room (Vacuumschmelze 
GmbH, Germany) using a 151 channel whole-head MEG 
system (CTF systems Inc., Canada) [19]. A third order 



Table 1 Characteristics of the included subjects 





CT + CRT 
(n = 14) 


CT 
(n = 18) 


Controls 
(n = 35) 


Gender, male % (n) 


57.1 (8) 


44.4 (8) 


45.7 (16) 


Age at assessment, y, M (SD) 


31.0 (4.3) 


24 (2.9) 


26.6 (6.1) 


Age at diagnosis, y, M (SD) 


5.7 (3.4) 


4.1 (2.0) 


N/A 


Time since diagnosis, y, M (SD) 


25.3 (2.5) 


20.0 (2.0) 


N/A 


Estimated IQ, M (SD) 


91.9 (18.1) 


105.8 (19.9) 


108.8 (19.4) 



IQ was estimated using a 4-subtest short-form of the Wechsler Adult 
Intelligence Scale Revised (WAIS-R III). M mean; SD standard deviation. 



software gradient was used with a recording passband of 
0.25-125 Hz and a sample frequency of 625 Hz in 49 
subjects (assessment around 2009), and a sample fre- 
quency of 312.5 Hz in 18 subjects (assessment around 
2007) [20]. Both sample frequencies occurred in each of 
the subject groups. Head position was monitored. MEG 
was recorded during a no-task, eyes-closed resting state 
condition. For this study, 149 of the 151 channels could 
be used. Additionally, 7 channels appeared to produce 
no signal for some of the subjects and were appropri- 
ately excluded from the calculation of the mean global 
and regional powers per subject. MEG recordings were 
converted to ASCII files and from these files four arte- 
fact free epochs of 13 s per subject were carefully 
selected by visual inspection by the two first authors 
(based on consensus). The MEG data were filtered in 
the following frequency bands: 6 (0.5-4 Hz), 6 (4-8 Hz), 
al (8-10 Hz), a2 (10-12 Hz), [3 (13-30 Hz) and y (30- 
50 Hz). Hypothetically, the lower and upper a bands are 
associated with different cognitive processes and many 
studies have found differential effects in these subbands 
[21,22]. There is no such rationale for subdividing the 
other frequency bands. Data in the 5 band were largely 
corrupted by movement artefacts like breathing. There- 
fore, all subjects with artefacts in the 8 band were 
excluded and group differences were investigated. 
Within this artefact-free dataset (CT + CRT n = 10, CT 
n=ll, controls n=12), the groups did not differ on 6 
power. Therefore, the 5 band was disregarded and more 
subjects with artefact-free data in the other frequency 
bands were included. Eventually, only 16 subjects (CT + 
CRT n = 2, CT n = 5, CON n = 9) were excluded due to 
movement artefacts in the other frequency bands, mostly 
caused by eye blinks and heart beats. 

Relative instead of absolute powers were analyzed, be- 
cause absolute powers are partially determined by the 
distance between the MEG sensor and the measured 
neural substrate. This would result in misleading out- 
comes, because part of our population has a smaller 
head size due to CRT [23]. The data were processed 
using CTF DataEditor and BrainWave 0.8.80 [19,24]. 
Relative band powers were computed for four 13 s 
epochs for each subject in the five frequency bands, and 
averaged for each subject. MEG channels were clustered 
based upon the approximate underlying cortical areas 
(see Figure 1). Regional means were calculated averaging 
relative powers of the available channels per cluster, over 
four epochs. 

Neurocognitive assessment 

The Amsterdam Neuropsychological Tasks (ANT) pro- 
gram was used to assess EF [25]. The computerized ANT 
provides for highly standardized assessments and auto- 
mated recordings of speed and accuracy of information 
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Figure 1 Schematic projection of MEG channels. Head drawing 
to illustrate the clustering of channels per region. Midline channels 
were disregarded. L = left, R = right, F = frontal, C = central, 
P = parietal, 0 = occipital, T = temporal. 



processing, attention processes and working memory. It 
has proven to be helpful in defining neurocognitive def- 
icit profiles in various clinical domains associated with a 
generally diffuse impact on the brain, such as phenylke- 
tonuria, multiple sclerosis, neurofibromatosis, and mid- 
dle late effects of childhood ALL [12,13,26-31]. Based on 
these studies, outcome parameters of baseline response 
speed, complex visual information processing, sustained 
attention, work pace and attentional fluctuations, cogni- 
tive flexibility (set shifting and inhibition), visuomotor 
skills, and visuospatial sequential working memory were 
selected for assessment of this study's population. Short 
task descriptions are given in Table 2, and in more detail 
in Additional file 1. The reliability and validity of these 
tasks have extensively been described elsewhere [32]. 
Task parameters that discriminated between the groups 
were selected for correlation with MEG data. 



was computed to estimate effect sizes (weak effect: n p 2 
«0.03; moderate: n p 2 «0.06; large: n p 2 >0.14) [34]. 
Cohen's d was used for simple contrasts (small effect: 
0.2 < d < 0.3; medium: 0.3 < d < 0.8; large: d > 0.8). 

Correlations between regional powers and discriminative 
cognitive variables 

The ANT variables that discriminated between patients 
and controls were correlated with MEG regional power 
outcomes in the discriminating frequency bands. Pear- 
son's correlations (r) and their p-values were calculated 
between regional powers and these neuropsychological 
variables, controlling for age. When r in the patient 
group or in the control group was significant, the r- 
values were compared. First, they were transformed into 
z'-values according to the Fisher transformation [35]. 
The difference |z*| between ri within the patient group 
and r 2 within controls was calculated and significance of 
that difference was determined by comparing |z*| to 
normal z-values [35]. 

Linear regression models of powers predicting cognitive 
performance 

Separately for each frequency band, regional means of 
relative powers that differed significantly between groups 
were entered into linear regression analysis predicting 
cognitive parameters that discriminated between groups, 
together with age. The backward method was used, 
using a probability (F) of 0.05 for entry and 0.10 for 
removal. 



Statistical analyses 

Individuals' global and regional means of relative power 
were entered in SPSS software (version 16.0; SPSS Inc, 
Chicago, IL) for statistical analyses. All powers were 
transformed according to the following logarithmic 
function [33]: 

( power \ 

power \n= log 

\1 — power] 

Global relative power differences between subject 
groups were analyzed using a Multivariate General Lin- 
ear Model (GLM) with group as between-subjects factor 
and the frequency bands as dependent variables. Age 
was used as a covariate when significantly correlated 
with the dependent variable within the healthy control 
group. Wilks' Lambda corrected F and p-values are 
reported. For each frequency band, regional power dif- 
ferences between groups were analyzed using a MAN- 
OVA for ten regions (bilateral central, frontal, parietal, 
occipital and temporal power). Differences between 
groups were tested with simple contrasts comparing 
each patient group with controls. Partial eta squared (r|p) 



Results 

Global relative power 

After e log transformation, separate inspection per group 
revealed normal distributions within all frequency bands, 
except for al in controls (SW statistic: 0.926 (p = 0.022)). 
Based on these results, the e log transformation was con- 
sidered sufficiently effective to allow for parametric stat- 
istical testing. 

The groups differed significantly in age (p = 0.001) due 
to the consecutive time periods the treatment protocols 
were applied. The correlation analyses within the control 
group between the powers and age resulted in a signifi- 
cant outcome for 9_ln (r = -0.560, p = 0.0005), so age was 
used as a covariate in the multivariate power analyses. 
The GLM results for the mean global relative powers are 
shown in Table 3. In the CT + CRT group, 0 power was 
slightly increased compared to controls and ct2 power 
was significantly decreased. Within the CT group, global 
power levels did not differ from controls. To facilitate in- 
terpretation, Figure 2 displays the untransformed mean 
global relative power values per group for all frequency 
bands. For additional graphs of transformed and untrans- 
formed global relative powers, see Additional file 2. 
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Table 2 Amsterdam Neuropsychological Tasks program: description of subtests 



TASK 



VARIABLE 



LABEL 



Baseline speed T_bs 
S_bs 

Feature Identification Ts_fi 
Td_fi 
diffLfi 

Memory Search Objects 2D T1_2d 
T2_2d 



Sustained attention 



Shifting attentional set 



Tracking 



Pursuit 



Visuospatial sequencing 



diffT_2d 

T_sa 

SD_sa 

PM_sa 

PF_sa 

TJnhib 

PJnhib 

T_flex 

PJlex 

Da_tr 

S_tr 

Dal(r)_tr 

D_pu 

S_pu 

Dl(r)_pu 

Sl(r)_pu 

Nit_vs 

Nitco_vs 



Baseline speed (simple reaction time) 
Baseline speed stability 

Speed of processing complex visuospatial information (similar condition) 
Speed of processing simple visuospatial information (dissimilar condition) 

Ts_fi minus Td_fi; Extra response time needed for higher order complex visuospatial information processing 

Speed of working memory search processes (low memory load condition); Target detection 

Speed of working memory search processes (high memory load condition); Requires continuous 
monitoring and updating of the contents of the working memory 

T2_2d minus Ti_2d; Extra response time needed for memory search processes 

Work pace; Response speed on sustained attention 

Attentional fluctuations; Response speed stability 

Percentage of misses (errors on target signals) 

Percentage of false alarms (errors on nontarget signals) 

Speed of inhibition 

Accuracy of inhibition (% errors) 

Speed of set shifting (flexibility) 

Accuracy of set shifting (flexibility)(% errors) 

Visuomotor accuracy; Mean absolute distance to ideal trajectory, mean left and right hand 
Visuomotor stability; mean left and right hand 
Visuomotor accuracy; left (right) hand 

Visuomotor accuracy; Distance to randomly moving target, mean left and right hand 
Visuomotor stability, mean left and right hand 
Visuomotor accuracy; left (right) hand 
Visuomotor stability; left (right) hand 

Visuospatial working memory; Number of correctly identified dots in a visuospatial pattern 
Sequential visuospatial working memory; Number of correctly identified dots in correct order 



diff_Nit_Nitco Nit_vs minus Nitco_vs; Accuracy of sequential working memory processes 



For more elaborate descriptions of these subtests, see Additional file 1. 

Regional relative power 

For each frequency band, a separate MANOVA was used 
to test for regional power differences between groups in 
the ten regions shown in Figure 1. The multivariate tests 
did not reach significance, with just a trend for the 

Table 3 Global relative power differences 



group difference on regional a2 powers (p = 0.067). 
Group contrasts versus controls within the 9 band 
revealed increased powers for the CT + CRT group in 
the RC (p = 0.038) and RP (p = 0.029) regions and trends 
for increases in the LO (p = 0.066) and RO (p = 0.051) 



MANOVA 



Between-Subjects effects 



Simple contrasts (versus controls) 





AGE 


Sign. 
(P) 


Effect 
size 
(Hp) 


GROUP 


Sign. 
(P) 


Effect 
size 
tip) 


CT + CRT group 
Sign, (p) Cohen's d 


CT group 
Sign, (p) Cohen's d 


Multivariate tests 


F[5,59] = 2.S16 


0.039 


0.176 


F[10,1 18] = 1.809 


0.066* 


0.133 


N/A 


N/A 


N/A 


N/A 


9_ln 


F[1 ,63] = 12.404 


0.001 


0.164 


F[2,63] = 2.191 


0.120 


0.065 


0.069* 


0.190 


0.582 


0.061 


aljn 


F[1,63] = 1.589 


0.212 


0.025 


F[2,63]= 1.057 


0.354 


0.032 


0.152 


0.335 


0.693 


0.201 


a2_ln 


F[1 ,63] = 4.404 


0.040 


0.065 


F[2,63] = 5.757 


0.005** 


0.155 


0.006** 


-0.719 


0.243 


0.207 


PJn 


F[1,63] = 1.509 


0.224 


0.023 


F[2,63] = 0.511 


0.602 


0.016 


0.889 


0.171 


0.347 


-0.371 


Y_ln 


F[1, 63] = 1.810 


0.183 


0.028 


F[2,63] = 1 .684 


0.194 


0.051 


0.117 


-0.386 


0.218 


-0.448 



Differences were tested between the CT + CRT group (n = 14), the CT group (n = 18} and controls (n = 35) using multivariate GLM, with age as a covariate. 
* = statistical trend, ** = significant difference. 
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Figure 2 Mean global relative power values. X-axis: Frequency band; y-axis: Mean global relative power (%). Error bars represent one standard 
deviation. 



regions. This group also displayed two trends for 
increased levels of al power in the LO (p = 0.081) and 
RO (p = 0.059) regions. Within the a2 band, the CT + 
CRT group displayed decreased power levels in all 
regions, except for LF. The CT group only displayed a 
significant increase of ct2 power (p = 0.030) in the left 
frontal region compared to controls. For a table of all 
statistics on the multivariate tests and contrasts, see 
Additional file 3. Significant contrast results for the CT + 
CRT group versus controls are displayed in Figure 3. 
The directions of the deviations are colour coded per re- 
gion (red for increases, blue for decreases). No signifi- 
cant regional Group x Hemisphere interactions were 
found. 

The time that has passed since treatment could hypo- 
thetically allow for either recovery or worsening of late 
effects. Therefore, time since diagnosis was correlated 
with powers in all frequency bands, separately within 
each patient group. There were no significant correla- 
tions with either global or regional powers. Hypothetic- 
ally, a younger age could mean more vulnerability to 




■ Power higher in the CT+CRT group at p < 0.05 
Power higher in the CT+CRT group at p < 0.10 

■ Power lower in the CT+CRT group at p <0.05 

■ Power lower in the CT+CRT group at p<0.10 

Figure 3 Visualization of regional power differences between 
the CT + CRT group and controls. Differences concern the 8 and 
a2 frequency bands. The directions of the deviations are colour 
coded per region (red for increases, blue for decreases). 



neurotoxicity of treatment. However, no effect of age at 
diagnosis was found on either global or regional powers. 

Group differences on cognition 

The cognitive variables are explained in Additional file 1. 
Age did not correlate significantly with any of the cog- 
nitive parameters within the control group and was 
therefore not used as a covariate. The CT + CRT group 
performed worse than controls on measures of cognitive 
flexibility (percentage of errors (P_flex, p = 0.002), atten- 
tional fluctuations during sustained attention (SD sa, 
p = 0.062), visuomotor accuracy (Dpu, p = 0.049; S_pu, 
p = 0.075; Dr_pu, p = 0.015; Sr_pu, p = 0.029) and se- 
quential visuospatial working memory (Nit_vs, p = 0.065; 
Nitco_vs, p = 0.023; Nit-Nitco_vs, p = 0.036). A table of 
statistics, including effect sizes, can be found in Add- 
itional file 4. When the CT group was compared with 
controls, trends were found for worse sequential visuo- 
spatial working memory (Nitco_vs, p = 0.092) and accur- 
acy of inhibition (P inhib, p = 0.059). 

Correlations between regional powers and discriminative 
cognitive variables 

Correlations were calculated within the composite group 
of CT + CRT patients and controls to investigate general 
patterns. Significantly differentiating regional powers 
were correlated with differentiating cognitive variables. 
Increased levels of RC and RP 0 power correlated signifi- 
cantly with worse performance on visuomotor control 
(D pu, Dr pu and Sr pu). Decreased oc2 powers in the 
eight significant regions of the CT + CRT group (LO, LP, 
LT, RC, RF, RO, RP and RT) were significantly correlated 
with better performance on attentional fluctuations 
(SD sa) and with worse performance on visuomotor 
control (D_pu and Dr_pu). Correlations within the CT + 
CRT group and within controls were separately calcu- 
lated and then compared (see Additional file 5). There 
was a trend for the difference between CT + CRT 
patients and controls in correlation coefficients of 
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attentional fluctuations with right frontal al powers. 
This correlation was absent in the CT + CRT group 
(r = 0.078), while it was significantly positive in controls 
(r = 0.549). Correlations between visuomotor accuracy of 
the right hand (Drpu) and regional al powers differed 
significantly in seven regions (not in RF) between CT + 
CRT patients (where there were significantly negative 
correlations with r ranging from -0.508 to -0.697)), and 
controls (where the correlations were not significant and 
ranging from -0.129 to 0.142). Roughly the same pattern 
was seen in the 8 band, but with positive correlations 
within the CT + CRT group between Dr pu and RC 6 
(r = 0.630) and RP 6 (r = 0.599). 

To verify whether the increased LF al power level in 
the CT group was associated with the mild deficiencies 
in accuracy of inhibition and sequential visuospatial 
working memory, these outcomes were correlated, but 
no significant results were found. The significant corre- 
lations between oc2 and attentional fluctuations observed 
in controls were absent in the CT group. 

To explore the possibility that increased levels of 0 
power reflect the results of compensatory activity, the 
data were searched for associations between increased 
regional 6 powers and better cognitive performance on 
any of the neuropsychological tasks. Two trends were 
found for increased left occipital 6 power in association 
with less attentional fluctuations (r = -0.509, p = 0.075) 
and better right hand visuomotor accuracy on the track- 
ing task (r = -0.509, p = 0.075). A description of the 
tracking task can be found in Additional file 1. 

Linear regression models 

The differentiating regional powers were entered in lin- 
ear regression analysis for each group separately, and for 
the composite group of CT + CRT patients and controls. 
Of the three pursuit variables, visuomotor accuracy of 
the preferred hand (Dr pu) is best explained by regional 
0 power. Within the CT + CRT group, 39.6% variance on 
Dr pu could be explained by increased RC 0 power. No 
significant model could be reached for controls. 

Attentional fluctuations (SD_sa) could be significantly 
explained by oc2 powers in one or two regions, regardless 
of group composition. More variance could be explained 
in the CT + CRT group (R 2 = 0.653), than in controls 
(R 2 = 0.345). D_pu and Dr_pu could be significantly 
explained by al powers in multiple regions. Similarly, 
more variance could be explained in CT + CRT patients 
(R 2 = 0.912 (D_pu) and 0.727 (Dr_pu)) than in controls 
(R 2 = 0.173 (D_pu) and 0.327 (Dr_pu)). An overview of 
linear regression outcomes is given in Additional file 6. 

Discussion 

The global power spectrum of CRT-treated patients 
reflects a slight increase of 0 activity and a significant 



decrease of al activity. Patients treated with CT alone 
display a power spectrum similar to controls, except for 
a significantly increased level of left frontal al power. 

Late effects of cranial irradiation 

The increase in 0 and decrease in al indicate a tendency 
towards global slowing of brain oscillatory activity, a 
finding that is consistent with other recent MEG studies 
in low grade glioma patients and Alzheimer patients 
[15,16]. Increased 0 power is also known as a nonspeci- 
fic sign of brain pathology [17]. Our results confirm 
Tucker's findings from 1989 of increased 0 and 
decreased a power in irradiated patients and are also in 
line with results from a review by Klimesch, published in 
1999 [18,36]. He concluded that a power is lowered and 
0 power enhanced in subjects with a variety of different 
neurological disorders. Also in healthy people in the late 
part of their lifespan this pattern is familiar and is sup- 
posed to reflect some degree of normal aging in the 
form of vascular or fibrillary degeneration [17]. Taken 
together with the fact that dementia is frequently 
reported as a late effect of whole brain irradiation and 
the neuropsychological deficiencies currently present, 
these findings suggest that the irradiated brain might be 
aging faster and could be at risk for early-onset dementia 
[18,37-40]. Long-term follow-up of this patient group is 
therefore increasingly important. Longitudinal research 
is also recommended. 

High percentages of variance of visuomotor accuracy, 
especially of the right hand, are explained by increased 0 
and decreased al powers in the CT + CRT group. This is 
a strong indication that these deviant power levels are 
underlying this particular deficiency. Unexpectedly, 
decreased regional al powers are significantly associated 
with less attentional fluctuations in both irradiated 
patients and controls. In general, a is suppressed by 
stimuli (e.g. light) and mental activities [17]. Rihs et al. 
(2009) describe how a power increases and decreases 
co-vary with visual cortex excitability related to anticipa- 
tory visual attention processes [41-43]. This illustrates 
that patterns of a power fluctuations are difficult to de- 
cipher and it is almost impossible to interpolate between 
resting state and attentional task data. 

The CT + CRT group demonstrates decreased al levels 
and worse attentional fluctuations. Within the right 
frontal region of irradiated patients, where the al power 
level is decreased, the association with attentional fluc- 
tuations is absent, while there is a strong positive associ- 
ation in controls. This suggests that especially the 
significantly decreased al power in the right frontal re- 
gion underlies the deficiency in sustained attention in 
the CT + CRT group. Lower al is not helping sustained 
attention here. The al associations in the occipital 
regions are slightly, although not significantly, stronger 
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than in controls, which could hypothetically indicate 
that the occipital lobes are trying to compensate for 
decreased a2 power in the right frontal region. 

Additionally, recent research on 9 activity, e.g. by 
Voytek et al. (2010) [44], suggests that increased 0 power 
may reflect compensatory activity for a damaged brain 
region that is being challenged, rather than reflecting 
pathology. In our study, subtle implications were found 
for compensatory 0 activity in the left occipital lobe. 
Combined, the increase in left occipital 0 power and 
the decrease in occipital oc2 powers both being asso- 
ciated with less attentional fluctuations, suggest that at- 
tentional fluctuations in irradiated patients may be 
partly counterbalanced by compensatory activity in the 
occipital lobes. One might arguably postulate that the 
observed trends for elevated occipital al power repre- 
sent nothing more than a side effect of working with 
relative powers instead of absolute. These slightly 
increased levels were not associated with any of the 
cognitive deficiencies in the irradiated patients, although 
reports have been made about associations between an 
increase in lower a power and disturbed working 
memory [15]. 

Late effects of chemotherapy 

The patient group treated with CT alone showed no signs 
of early aging. The subtlety of the cognitive deficiencies 
and limited spectral power changes are consistent with 
overall findings that the patients treated with CT only do 
not perform significantly worse than controls. The sig- 
nificant reduction of late effects makes CT a superior al- 
ternative to CRT, especially since there was no trade-off 
in terms of mortality or recurrence [1,2]. Event-free sur- 
vival at 10 years after treatment with ALL-6 was even 
30% better than after treatment with ALL- 5 and CNS re- 
lapse rates improved from 12.9% to 1.1% [45]. Mean- 
while, suspected neurotoxic effects of CT might be 
masked by compensatory mechanisms. It remains un- 
clear whether the elevated level of left frontal a2 power 
should be interpreted as compensatory or pathological. If 
lowered a2 is regarded as compensatory, increased a2 in 
the CT patients should be regarded as pathological. 
However, increased a2 power and the mild cognitive de- 
ficiencies observed in the CT patients were not asso- 
ciated, which contradicts the pathology hypothesis. In 
CT patients, oc2 power is not correlated with attentional 
fluctuations, so the absence of attentional fluctuations 
combined with increased oc2 is not necessarily inconsist- 
ent with the hypothesis of lower a2 levels being compen- 
satory. Also, the possibility of increased oc2 being 
compensatory in other cognitive domains cannot be 
ruled out. 

The time that passed since end of treatment (20- 
25 years) and the age of the patients at time of treatment 



(an age of relatively high plasticity of the brain) could 
have allowed for alternative brain strategies to develop. 
For the CT patients, this could have resulted in mostly 
normal cognitive performance. For the irradiated group, 
however, neurotoxicity is much more severe. Compensa- 
tion might exist, yet appears to be inadequate, which 
emphasizes this group's need for long-term follow-up 
and interventions. 

Methodological considerations 

Although there were good reasons for analyzing relative 
powers instead of absolute, this also creates a limitation. 
A decrease in one frequency band might cause an in- 
crease in another band, although this would not be the 
observed pattern in absolute measures. On the other 
hand, consistency with the literature and significant 
associations with cognitive deficiencies support our find- 
ings in the CT + CRT group. Another possible limitation 
of the study is the relatively small sample size which 
raises the question whether this set of subjects is repre- 
sentative for the whole population. Confirmation by fu- 
ture studies is needed. 

No corrections were applied for multiple compari- 
sons, although commonly the Bonferroni method would 
be applied in order to lower the chances of false discov- 
ery of significant results (the type I error rate). In a 
publication by Perneger et al. (1998), a view widely held 
by epidemiologists is expressed that this method creates 
more problems than it solves and the best way of deal- 
ing with multiple comparisons is by simply describing 
what tests of significance have been performed, and 
why [46]. This way, findings are not interpreted differ- 
ently according to how many other tests were per- 
formed. Also, type I errors cannot be prevented without 
increasing the number of type II errors, which are 
equally worth avoiding. In the present study, pre- 
established hypotheses were tested in our correlation 
analyses, instead of just correlating everything, and more 
than 50% of these hypotheses were confirmed. In 
addition, the correlation results were supported by the 
regression results. This, combined with reporting effect 
sizes in addition to significance levels, should provide 
the reader with sufficient information to base an inter- 
pretation of the results on. 

Conclusions 

Power spectrum analysis of MEG registrations is able to 
demonstrate abnormal patterns of resting-state oscilla- 
tory brain activity in irradiated long-term survivors of 
ALL. The suspicion of a faster aging brain after CRT 
warrants careful long-term follow-up and screening for 
early-onset dementia. In future ALL therapy, CRT 
should be completely abolished. 
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Additional file 1: ANTprogram. Description of subtests of the 
Amsterdam Neuropsychological Tasks program (ANT). 

Additional file 2: Global relative powers. Mean global relative powers 
per group. Four visualisations of the same data. A: Bar chart of the raw 
power values, equal to Figure 2 in the main text. B: Bar chart of the log 
transformed power values. C: Line chart of the raw power values. D: Line 
chart of the log transformed power values. Error bars represent one 
standard deviation. 

Additional file 3: Regional relative power differences. Statistics of 
the multivariate tests and group contrasts of the regional relative powers 
in each frequency band. (MANOVA 1: theta, MANOVA 2: alpha 1, 
MANOVA 3: alpha 2, MANOVA 4: beta, MANOVA 5: gamma). Group 
differences were tested with age as a covariate. 

Additional file 4: Neuropsychological outcome. Overview of group 
differences on all neuropsychological variables assessed with the 
Amsterdam Neuropsychological Tasks (ANT) program. Separate ANOVA's 
were applied per variable, with simple contrasts between each patient 
group and controls. 

Additional file 5: Correlations. Overview of differences between 
correlation coefficients of regional power and ANT variables, 
calculated separately within the CT + CRT group and in controls. 

Additional file 6: Linear regression models. Overview of linear 
regression models within the separate (composite) groups. R 2 indicates 
the explained variance of ANT variables by regional powers; p is the 
significance of each linear regression model. 
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